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ABSTRACT 


A computer model based on the method of images was used to investigate and 
calculate the pressure amplitude distribution in a wedge-shaped (tapered) water layer 
overlying a fast, absorbing bottom. The pressure amplitude of the field was generated 
and compared at reciprocal positions of source and receiver. The results showed that the 
ratio of pressure-difference to the pressure when the position of source or receiver was 
more near the bottom, that is, deeper and more near the apex. 
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I. INTRODUCTION 


Propagation of sound radiated from a source in an ocean with a sloping 
bottom has received considerable attention during recent years for both theoretical 
reasons and applications to underwater sound. Many scientists [Ref. 1-11] have 
researched this subject and created a number of acoustic models to predict the 
sound field within a wedge-shaped fluid overlying a penetrable or a rigid fluid 
bottom. 

Considering practical interest for underwater acoustics, a limitation of the 
ideal wedge model of the ocean lies in the boundary conditions: The top boundary 
may be approximated as an impenetrable, pressure-release surface, but the bottom 
of a real ocean wedge is not usually a perfect reflector. It can be considered as 
a fast bottom, that is, a water fluid sediment interface that shows a critical grazing 
angle. 

The acoustic field in a wedge shaped shallow water duct with ideal boundary 
conditions has been studied by Bradley and Hudimac [Ref. 1]. They have analyzed 
the case of an isospeed duct with one pressure release and one rigid surface. The 
theoretical analysis has been carried out in both image theory and normal mode 
theory. 
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pressure amplitude and phase in the upslope direction along the bottom of wedge- 
shaped fluid layer overlying a fast fluid bottom. 

Jensen and Kuperman [Ref. 4] used the parabolic equation model to study the 
model cut off during upslope propagation in a wedge-shaped ocean and made 
comparison between e3q)eriment and theory in 1980. 

In 1982, Lee and Botseas [Ref. 5] used an Implicit Finite-Difference (IFD) 
computer model, that is problems arise when the fourier transform encounters an 
interface between two medium having different sound speed and densities, to solve 
the parabolic equation. 

In 1983, Jager [Ref. 6] developed IFD program that incorporates exact 
interface conditions which are preserved for horizontal and sloping interfaces along 
a user-specified bottom profile for solving the parabolic equation. 

The sound field in an absorbing fluid substrate underlying a wedge-shaped 
fluid with high sound speed has been studied by Coppens, Sanders, and Humphries 
in 1984. [Ref. 7] 

In the same year, Baek [Ref. 8] used the computer model based on the method 
of image to predict pressure amplitude and phase distribution in a wedge-shaped 
medium overlying a fast absorbing bottom. 

Lesesne [Ref. 9] also studied and compared two computer models developed 
by Coppens and Sanders which is not limited to up or downslope direction. 
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In this research, I will use computer model for image method that was 
programmed in Fortran for use on the NFS IBM 3033 computer main frame to 
calculate the pressure field in a wedge-shaped ocean with a pressure release 
surface and an acoustically fast bottom. 

This computer program was developed by LTJG George Nassopoulos(Greece 
Navy). The ou^uts will be compared at reciprocal positions of source and 
receiver and be analyzed for the evaluation of the validity of the image model in 
above conditions. 
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n. THEORY 


The following development is a summary of previous research of J.V.Sanders, 
A.B.Coppens and their students at the Naval Postgraduate School. 

The complex acoustic pressure in the wedge from the point source can be 
determined by using the method of images. The assumption of constant sound 
speed in each of the two different mediums makes the method of images an 
appropriate approach for predicting the sound field. 

In the method of images, the reflections of sound from boundaries are replaced 
by images of the source that reproduce the effects of the interactions with the 
boundaries. For the wedge-shaped layer, images He on a circle whose center is the 
apex of the wedge . The source and each of image radiates a spherical wave of the 
appropriate phase. The total pressure and phase at any field point within the wedge 
is found by the phase-coherent summation of these waves. Higher images 
correspond to more reflections from the bottom so the pressure distribution of 
these higher images decrease progressively . 

A typical three-dimension geometry for wedge is shown in Figure 1. The 
foUowing definition and symbols for source and receiver parameters apply and wiU 
be used throughout. 
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Ri 

R2 

Yo 

7 

5 

D1 

CC 

e. 

X 

(X/K.2 = 


wedge angle 

normalized distance from apex to source 

normalized distance from apex to receiver 

normalized distance along the shore between receiver and 

source 

source angle measured upward from the interface 

receiver angle measured upward from the interface 

the ratio of medium density to bottom density (P 1 /P 2 ) 

the ratio of the speed of sound in the medium to the speed of 

sound in the bottom (C 1 /C 2 ) 

critical angle for the bottom 

normalizing distance 

wave number in the bottom divided into the absorption in the 
bottom (loss term) 


For the fast bottom, all distances are normalized to the dump distance x 
measured from the apex along the wedge interface at which the lowest mode 
attains cut off. The following formula is applied, 
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Figure 1. Three-Dimensional Wedge Geometry 


K,X= -—- 

^ 2sm0.tanP 


The normalized range is used in the image method program. For this model, 
we use wedge angles of the from x/N where N is an integer and consequently no 
apex term exists. [Ref. 2] 

Then, the number of images in upper or lower half-space is 







In Figure 2, the relationship between the receiver and the n*"* source is 
represented and the angle of the n*** image, measured from the bottom, is given 
by: 

0^=P(n-l)+Y , for n odd 
» for n even 

Figure 3 shows the distance between the receiver and the n* image; 
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(R'f^R^^^^^R2^~2R^R2Cos(Q^-b} 


R(n)=^l(Rf^=^lRf^ 

for the upper group of images, and 

(Rf=R^^*R2^-2R^RJCOsiQ^+b) 


RXn)=^R^^+R^^-2R^RjCos(d^+b)+YQ^ 
for the lower group of images. 
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The number of times of the relevant reflection along the path from image to 
receiver is equal to the number of times the path interacts the bottom and its 
images. The total reflection coefficient is the product of all reflection coefficients 
along the path from the specific image to the receiver. 

The reflection coefficient for the nj^ image at the m^, bounce can be 
approximated by the plane wave reflection coefficient, if the distance from the 
source to the point on the bottom at which the initial bounce takes place is large 
compared to wave length in the wedge. 

The reflection coefficient is; 




92^2 

Pi^i 






nm 


nm 


where 






i-(—)cos2e 


nm 


sin6 


nm 
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With a lossy bottom, the above ean be generalized with the help of 


^1 ^1 tL 

^2 C^a-A ^2 A 
K 


and 


/* 2 /* 2- /« 2 

A <-) *2jA f 

^2 ^2 ^2 ^ 


The result[Ref. 15] is 


^nm~ 


h 

Pi 


sin6 

^ ^nm 


(^)We. 


—sin0^+ 
Pi 


(—) -cos^O 

N ^2 


nm 


—sm0 — +j—^)f^^-a 
Pi _ 

—sin0/im+—/^^i> 2 +^ 2 ~ 
Pi 


where 


«.(^) -cos^e„. f 

^2 ^2 ^ 


So, the angle of incidence must be determined to find a reflection coefficient. 
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The angle of incidence of the upper n”* image on the m*^ plane is shown in 
Figure. 4: 




/?jSiii(0^-2i» P) -'-/{2Sin(2m P - 6) 

m 


m=l, 2 , 


for the upper image, 


SHeJ= 


i{^sin(0^-2mp)+/^in(2mp +6) 


m=0y 1, 2, ...M 


for the lower image. 

In these expressions the index m indicates each bottom plane and the number 
M is the total number of bottom the sound interact with along path from source 
to receiver; 


M=INT [-? 5 ] = INT [—] 

P 2 

The complex acoustic pressure at the arbitrary position is the summation of the 
complex pressure from every image which consist of the contribution from the 
direct path, the images above the wedge and the image below the wedge. In this 
program, upper family of images include the direct path. So, we don’t need to 
calculate the reflection coefficients for the first and second images among upper 
images because these don’t interact with any bottom. 
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Figure 4. The Angle of Incidence 6 ^^ of the Upper n^ Image 


on the m® Plane 


The complex pressure from the upper images result in; 


»=1 JiW B,=l 


Where 


= -1 for n - 2 


and from the lower images 




^1 

p.=y—L-exp 

R'in) 


n ^(0 

m»0 


The total complex acoustic pressure P is the sum of the P„ and P, 


m. EVALUATION AND DISCUSSION 


The effects of changing the source and receiver position on the prediction of 
the pressure field were evaluated at each reciprocal position. 

The following fixed parameters were used, 


0 

0 

0 

D1 

0.9 

CC 

0.9 

oj/Kj = 

0.0183 


The depictions of the pressure fields for these conditions as functions of Ri, 
R 2 , 7 , 8 and are shown in figures. The specific data are shown in tables. The 
model was tested by varying either source or receiver and holding the other 
constant in each Yo. The test was repeated after exchanging the positions of 
source and receiver. All the tests were done within ten normalized distances from 
apex. 

A. ANGLE OF SOURCE, RECEIVER 

The first case was examined by varying the angle with the source and receiver 
on axis-positions(Yo =0). As depicted in Figure 5, there are two positions, one 
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for source and another for receiver, and vise versa. The first position is set as 


constant R=10. The second position is set to R=1 at mid-depth 



The depiction of the pressure field for these conditions as function of angle is 
shown in Figure 6. The corresponding overplot is shown in Figure 7, and Table 
1 provides data of overall pressure. 

Comparisons of pressure amplitudes at each reciprocal position shows that 
there are 1 %-2.5% difference near the bottom (Figure 9 and 10), while they agree 
well near the surface (Figure 8) and at the middle depth. As the depth is deeper 
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near the bottom,the ratio of difference to pressure (AP/P) is larger. This pressure 
difference is probably the result of using the plane wave Rayleigh reflection 
coefficients, discussed by Kinsler, Frey, Coppens and SandersfRef. 10], when is 
only a good approximation if the sound source is not too close to the bottom. 

B. DISTANCE FROM APEX 

In this case, the distance from apex was varied with =0 and the same 
angular plane. One position is set as constant R=10 and another is varying from 
apex to R=l. The consistency of the pressure amplitude at each reciprocal 
position for various distance from apex is good near the surface and the mid-depth. 
But there are some pressure differences near the bottom. The overplot and data 
for this case are provided by Figure 11 and 12, Table 4. 

As the distance from apex decreased, the ratio of difference to pressure (AP/P) 
increased. 

C. SHORE DISTANCE AS A VARIABLE 

Until now, all studies were for on axis position (Yo—0). In the following 
discussion the shore distance (Yo) is the variable. All the above cases were re¬ 
examined with Y^ =5 and Y^ =10. 
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The graphs shown in Figure 13 through Figure 18 indicate the results of 
Yo =5 with a solid line and Yo =10 with the symbol connected by dashes. All the 
specific data were provided by Table 5 through 7. The pressure difference (AP) 
and the ratio of difference to pressure (AP/P) near the surface for various angles 
were shown in Figures 13 and 14. 

The results (Figures 15 and 16) show the ratio of difference to pressure (AP/P) 
still shows decreasing error as source or receiver move away the bottom and 
increases with increasing Y^,. 

Figures 17 and 18 show that as either the source or receiver approaches the 
apex, the ratio of difference to pressure(AP/P) becomes very high. 


17 


IV. CONCLUSION 


A computer program used was created to predict the pressure of a point source 
everywhere within a wedge-shaped fluid overlying a fast bottom using the image 
method. This thesis was an attempt to compare the pressure amplitude for various 
reciprocal positions of source and receiver and to evaluate the validity of the image 
model in these conditions. 

Figure 19 and 20 show the overplot of results for various configuration. Each 
graph indicates the overplot of pressure at reciprocal positions between source and 
receiver. =0 is indicated by a solid line, Yo =5 with the dashes and =10 
with the symbol connected by dashes. 

Comparisons of pressure amplitude within the wedge at Yo=0, 5, 10 for 
various angles are shown in Figure 19. Comparisons of pressure amplitude at mid¬ 
depth for various distances from the apex are shown in Figure 20. These 
comparisons show good consistency in each case. 

The pressure amplitudes are more consistent at reciprocal positions when the 
position is near the surface and at mid-depth. The ratio of difference to pressure 
increases as the bottom is approached. This difference is probably caused by the 
failure of the plan-wave reflection coefficient when the source or receiver is near 
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the bottom. 


In the case of near the bottom, the distance from apex affects the consistency 
of the pressure, except for near the surface and mid-depth cases. As the distance 
from the apex increases, the ratio of difference to pressure (AP/P) decreases. The 
ratio of difference to pressure (AP/P) is also proportional to the shore distance in 
near the bottom case. 

From the above results, it is concluded that the image model is valid for 
predicting the pressure in the sound field in a wedge-shaped layer overlying a fast 
penetrable bottom, except very close to the bottom and apex. 
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APPENDIX A: TABLES 


Table 1. PRESSIJRE WITfflN WEDGE FOR VARIOUS ANGLES (Y# = 0) 


7 

P 

Constant 

6 

P 

Relative 

0.0 

0.04913 

Rj = 10 

0.0 

0.04795 

Ri = 1 

1.0 

0.16694 

Ra = 1 

1.0 

0.16561 


2.0 

0.25990 

09 

II 

o 

2.0 

0.25882 

7 = 5V 

3.0 

0.31425 

o 

II 

3.0 

0.31365 

Yo = 0 

4.0 

0.32956 


4.0 

0.32930 


5.0 

0.31508 


5.0 

0.31501 


6.0 

0.27836 


6.0 

0.27841 


7.0 

0.22428 


7.0 

0.22442 


8.0 

0.15702 


8.0 

0.15724 


9.0 

0.08078 


9.0 

0.08108 


10.0 

0.00000 


10.0 

0.00050 



Table 2. PRESSURE NEAR SURFACE FOR VARIOUS ANGLES (¥„ = 0) 


7 

P 

Constant 

8 

P 

Relative 

9.0 

0.08078 

o 

II 

Pi 

9.0 

0.08109 

R , = 1 

9.1 

0.07283 

II 

9.1 

0.07315 


9.2 

0.06484 

o» 

II 

o 

9.2 

0.06517 

7 = 5 ° 

0.3 

0.05682 

o 

II 

9.3 

0.05715 

Yo = 0 

9.4 

0.04876 


9.4 

0.04910 


9.5 

0.04068 


9.5 

0.04103 


9.6 

0.03257 


9.6 

0.03293 


9.7 

0.02444 


9.7 

0.02481 


9.8 

0.01630 


9.8 

0.01668 


9.9 

0.00815 


9.9 

0.00854 


10.0 

0.00000 


10.0 

0.00050 
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Table 3. PRESSURE NEAR BOTTOM FOR VARIOUS ANGLES (Y„ = 0) 


H 

P 

Constant 

6 

P 

Relative 

DEG 

AP 

AP/P 

0.0 

0.04913 

II 

o 

0.0 

0.04795 

Ri = l 

0.0 

0.00118 

0.02402 

0.1 

0.06151 

Ra = 1 

0.1 

0.06029 

R 2 ~ 10 

0.1 

0.00122 

0.01983 

0.2 

0.07381 

8 =5^^ 

0.2 

0.07256 

7 =5“ 

0.2 

0.00125 

0.01693 

0.3 

0.08601 

o 

II 

0.3 

0.08474 

Yo = 0 

0.3 

0.00127 

0.01477 

0.4 

0.09809 


0.4 

0.09680 


0.4 

0.00129 

0.01315 

0.5 

0.11002 


0.5 

0.10872 


0.5 

0.00130 

0.01182 

0.6 

0.12180 


0.6 

0.12049 


0.6 

0.00131 

0.01076 

0.7 

0.13340 


0.7 

0.13207 


0.7 

0.00133 

0.00997 

0.8 

0.14480 


0.8 

0.14347 


0.8 

0.00133 

0.00919 

0.9 

0.15598 


0.9 

0.15465 


0.9 

0.00133 

0.00853 

1.0 

0.16694 


1.0 

0.16561 


1.0 

0.00133 

0.00797 


Table 4. PRESSURE NEAR BOTTOM FOR VARIOUS DISTANCES (Y, = 0) 


Ri 

P 

Constant 

R2 

P 

Relative 

R 

AP 

AP/P 

0.0 

0.00000 

7 = 0.2“ 

0.0 

0.00003 


0.0 

0.00003 

1.00000 

0.1 

0.00067 

6 = 0.2“ 

0.1 

0.00070 


0.1 

0.00003 

0.04286 

0.2 

0.00238 

Ra = 10 

0.2 

0.00278 


0.2 

0.00040 

0.14388 

0.3 

0.00815 

0 

II 

0.3 

0.00891 


0.3 

0.00076 

0.08530 

0.4 

0.01967 


0.4 

0.02081 


0.4 

0.00114 

0.05478 

0.5 

0.03616 


0.5 

0.03761 


0.5 

0.00145 

0.03855 

0.6 

0.05328 


0.6 

0.05486 


0.6 

0.00158 

0.02880 

0.7 

0.06713 


0.7 

0.06867 


0.7 

0.00154 

0.02243 

0.8 

0.07670 


0.8 

0.07813 


0.8 

0.00143 

0.01830 

0.9 

0.08246 


0.9 

0.08373 


0.9 

0.00127 

0.01517 

1.0 

0.08548 


1.0 

0.08659 


1.0 

0.00111 

0.01282 
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Table 5. PRESSURE DIFFERENCE, RATIO OF DIFFERENCE TO PRESSURE NEAR 
_SURFACE FOR VARIOUS ANGLES (Y, = S, Y, = 10) 


DEG 

AP 

AP/P 

Constant 

AP 

AP/P 

Relative 

9.0 

0.00042 

0.00552 

Ri.2=1,10 

0.00067 

0.00875 

Ri,2=1,10 

9.1 

0.00043 

0.00626 

o 

II 

0.00070 

0.01011 

It 

o 

9.2 

0.00Q43 

0.00703 

6 = 5“ 

0.00074 

0.01197 

6 = 5“ 

9.3 

0.00044 

0.00821 

II 

0.00078 

0.01435 

o 

II 

IQQI 

0.00044 

0.00955 


0.00081 

0.01729 


9.5 

0.00045 

0.01169 


0.00084 

0.02138 


9.6 

0.00046 

0.01489 


0.00088 

0.02780 


mm 

0.00047 

0.02116 


0.00091 

0.03789 


9.8 

0.00048 

0.03055 


0.00094 

0.05735 


9.9 

0.00050 

0.06158 


0.00097 

0.11162 


10.0 

0.00075 

1.00000 


0.00105 

1.00000 



Table 6. PRESSURE DIFFERENCE, RATIO OF DIFFERENCE TO PRESSURE NEAR 
BOTTOM FOR VARIOUS ANGLES (Y, =5, Y, =10) 



DEG 

AP 

AP/P 

0.0 

0.00087 

0.02148 

0.1 

0.00089 

0.01780 

0.2 

0.00090 

0.01515 

0.3 

0.00092 

0.01338 

0.4 

0.00093 

0.01191 

0.5 

0.00095 

0.01088 

0.6 

0.00097 

0.01006 

0.7 

0.00098 

0.00929 

0.8 

0.00099 

0.00866 

0.9 

0.00100 

0.00812 

1.0 

0.00102 

0.00774 


Constant 


R,.2=1,10 
7 = 5° 

6 = 5*’ 

Yo = 5 


AP 

AP/P 

0.00071 

0.02536 

0.00072 

0.02081 

0.00074 

0.01797 

0.00076 

0.01592 

0.00077 

0.01419 

0.00080 

0.01317 

0.00081 

0.01206 

0.00082 

0.01115 

0.00084 

0.01052 

0.00085 

0.00987 

0.00087 

0.00942 



I 
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Table 7, PRESSURE DIFFERENCE, RATIO OVER DIFFERENCE TWO PRESSURE 
NEAR BOTTOM FOR VARIOUS DISTANCES (Y« = 5, ¥„ = 10) 


R 

AP 

AP/P 

Constant 

AP 

AP/P 

Relative 

0.0 

0.00002 

1.00000 

R,.2=1,10 

0.00004 

1.00000 

Ri,2=1,10 

0.1 

0.00000 

0.00000 

7 = 0.2* 

0.00013 

0.08844 

o 

d 

II 

0.2 

0.00013 

0.05603 

6 = 0.2* 

0.00023 

0.07395 

d = 0.2“ 

0.3 

0.00043 

0.07963 

II 

0.00036 

0.06679 

II 

O 

0.4 

0.00077 

0.06190 


0.00055 

0.06180 


0.5 

0.00105 

0.04370 


0.00081 

0.05625 


0.6 

0.00121 

0.03176 


0.00105 

0.04719 


0.7 

0.00123 

0.02404 


0.00119 

0.03781 


0.8 

0.00117 

0.01907 


0.00118 

0.02932 


0.9 

0.00108 

0.01577 


0.00109 

0.02316 


1.0 

0.00097 

0.01328 


0.00097 

0.01868 
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Figure 7. Overplot of Fig.6(a) and 6(b) 
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Figure 8. Overplot of Pressure in the Reciprocal Position Near the Surface at 
Y,=0 for Various Angles 
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Figure 13. Pressure difference(AP) Near the Surface at Y,=5,10 for Various 
angles 
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DISTANCE FROM APEX(R1.R2) 
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Figure 19. Overplot of Pressure Amplitude in the Reciprocal Position 
at Y^*p, 5,10 for Various angles 
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APPENDIX G: COMPUTER PROGRAM 


♦POPIMA 

* IMAGE THEORY-CROSS SLOPE POINT FIELD PRESSURE 

i)caic>{(>{ei{e4(}|c9ic}fe})ci|(%i|ci)(i(c3(c3)ci(Kitc>lcsk)(ii)(iltifci(c:|(>|t4<* 

INTEGER LILN,S1,S2,N1,K 

REAL*4 B,CC,C2,D,D1,D2,G,PI,P1,P2,Q1,R1,R2,T, 

* T4,T6,W0,W1,Y0,Y1,Y2,Z1,Z2,Z3,Z4,Z5,Z6, 

* T1(900),R8(900),R9(900),S(900),C(900),E(900), 

* F(900),Y,Z,R3,AL,PZ(900) 

REAL TQQ,TQQ1,TQQ2,TQQ3 
PI =ACOS(-1.0D00) 

WRITE(**)’INPUTB,G,D1,CC,R1,R2,AL,Y0’ 

READ(6,*)B,G,D1,CC,R1,R2,AL,Y0 

C CALL EXCMS{’FILEDEF 11 CLEAR’) 

DATA NOO/11/ 

DATANOU/10/ 

CALL EXCMS(’FILEDEF 11 DISK POP6 PLOTTER Al’) 

CALL EXCMS(’FILEDEF 10 DISK POP6 DATA Al’) 

D=-0.1 

* 

K=0 

DO 20 M =0,10 
D=D+0.1 
K=K+1 

WRITE(NOU,2) B,G,D,D1 ,CC,R1 ,R2,AL, YO 
2 FORMAT(3X,’B = ’,F4.1,/,3X,’G =’,F4.1,/,3X,’D = ’,F4.1,/,3X,’D1 = 

* %F4.1,/,3X,’CC = ’,F4.1,/,3X,’R1 =’,F4.1,/,3X,’R2 = ’,F4.1,/,3X 
♦,’AL = ’,F6.3 ,/,3X,’YO = ’,F4.1) 


C INPUT PARAMETERS 
C B = WEDGE ANGLE (DEG) = 180/Nl 
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C G = SOURCE ANGLE (DEG) 

C D = RECEIVER ANGLE (DEG) 

C Nl= # OF IMAGE POINTS 

C R1 = SOURCE DISTANCE (IN DUMP DISTANCES) FROM APEX 
C R2 = RECEIVER DISTANCE (IN DUMP DISTANCES) FROM APEX 
C YO = DISTANCE (IN DUMP DISTANCES) ALONG APEX 

C D1 = RH01/RH02 

C CC = C1/C2 

C AL = ALPHA/K2 

C A =# OF RECEIVER POSITIONS 
C D = 5.0 

C ]----.] 

C ICHOOSE SLOW OR FAST BOTTOM BY VALUE OF SPEED RATIO ] 
C ]CC = C1/C2 . ] 

C MAIN PROGRAM 

N1 = INT(180./B) 

T6 = 180./PI 
B = B/T6 
G =G/T6 
D = D/T6 
C2 = CC**2 
TQQ = TAN(B) 

C -----— 

C DECISION ABOUT SLOW OR FAST BOTTOM 
C T4 = PI/(2*SIN(ACOS(CC))’T'AN(B)) FOR FAST BOTTOM 
C T4 = PI/(2*TAN(ACOS(l/CC))*TAN(B)) FOR SLOW BOTTOM 
C---—---- 

c 

IF (CC.LT.l) THEN 
TQQl =ACOS(CC) 

TQQ2 = SIN(T(3Q1) 

ELSE 

TQQl = ACOS(l/CC) 

TQQ2 = TAN(TQQ1) 
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ENDIF 

TQQ3 = 2.’T'QQ2rrQQ 

T4 = PI/TQQ3 

Q1 = 1/DSQRT(2.0D00) 

D2 = (Y0*Y0)+(R1*R1)+(R2*R2) 

R3 = 2.*R1*R2 
SI = 1 

C I---I 

C I TfflS DO LOOP CALCULATES THE THETA(N) | 

C I AND THE IMAGE SLANT RANGES R8(N) AND R9(N) | 

DO 30 N = 1, N1 
IF(Sl.GT.O) T1(N) = (N-1)*B+G 
IF(SLLT.O) T1(N) = N*B-G 

51 = - SI 

R8(N) = SQRT(D2-R3*C0S(T1(N)-D)) 

R9(N) = SQRT(D2-R3*C0S(T1(N)+D)) 

30 CONTINUE 

C I SUM THE PRESSURE OVER ALL IMAGES I 
C I--- 

PI = 0.0 
P2 = 0.0 

DO 40 N = 1, N1 

52 = (-l)**(INT(N/2)) 

C 1----— -^-1 

C I REFLECTION COEFFICIENTS ALONG NTH UPPER PATH I 

W1 = 2*C2*AL 
II = INT((N-l)/2) 

DO50 1=1,11 

SO) = ABS(R1*SIN(T1(N)-2*I*B) 

& +R2*SIN(2*I*B-D))/R8(N) 

IF(S(I).GE.1)S(I)=1 
C(I) = SQRT(1-(S(I)*S(I))) 

T = S(I)/D1 

WO = (-C2+(C(I)*C(I))) 
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Y = SQRT((W0*W0)+(W1*W1)) 

Z =ABS(WO) 

IF(Y.LE.Z) Y = Z 
Y1 = Q1*SQRT(Y+W0) 

Y2 = -Q1*SQRT(Y-W0) 

Z1 = T-Y2 
Z2 = -Y1 

Z3 = Z1/(Z1*Z1+Z2*Z2) 

Z4 =-Z2/(Z1*Z1+Z2*Z2) 

Z1 = T+Y2 
Z2= Y1 

Z5 = Z1*Z3-Z2*Z4 
Z6 = Z1*Z4+Z2*Z3 
EG) = Z5 
F® = Z6 
50 CONTINtrE 

C I--—-_| 

C I PRODUCT OF REFLECTION COEFFICIENTS ALONG NTH UPPER PATH | 

C I-----I 

Z1 = 0 
Z2 = 0 
Z3 = 0 
Z4 = 0 
Z5 = 1 
Z6 = 0 

IF(N.LE.2.00) GOTO 110 
DO 60 I = 1,11 
Z1 = E® 

Z2 = F® 

Z3 = Z5 
Z4 = Z6 

Z5 = Z1*Z3-Z2*Z4 
Z6 = Z1*Z4+Z2*Z3 
60 CONTINUE 

110 Z1 = Z5 

Z2 = Z6 
T =T4*R8(N) 
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Z3 = COS(T) 

Z4 = -SIN(T) 

Z5 = Z1*Z3-Z2*Z4 
Z6 = Z1*Z4+Z2'»'Z3 
PI = P1+S2*Z5/R8(N) 

P2 = P2+S2*Z6/R8(N) 

WRrrE(6,*) ’I=M, ’ P1=’,P1,’ P2=%P2 
11 = 11+1 

I--_1 

I REFLECTION COEFFICIENTS ALONG NTH LOWER PATH | 


DO 70 I = 1, II 

S(I) = ABS(Rl*SIN(Tl(N)-2*0:-l)*B) 

+ R2*SIN(2*a-l)*B+D))/R9(N) 
IF(S(I).GT.l) S(I)=1 
C(I) = SQRT(1.001-S(I)*S(I)) 

T = S(I)/D1 

WO = -C2+C(I)*C(I) 

Y = SQRT((W0*W0)+(W1*W1)) 

Z = ABS(W0) 

IF(Y.LE.Z) Y = Z 
Y1 = Q1*SQRT(Y+W0) 

Y2 = -Q1*SQRT(Y-W0) 

Z1 = T-Y2 
Z2 = -Y1 

Z3 = Z1/(Z1*Z1+Z2*Z2) 

Z4 = - Z2/(Z1*Z1+Z2*Z2) 

Z1 = T+Y2 
Z2 = Y1 

Z5 = Z1*Z3-Z2*Z4 
Z6 = Z1*Z4+Z2*Z3 
E(I)=Z5 
F(I)=Z6 
CONTINUE 


c 

c 

c 


i PRODUCT OF REFLECTION COEFFICIENTS ALONG NTH LOWER PATHj 


Z1 = 0 
Z2 = 0 
Z3 = 0 
Z4 = 0 
Z5 = 1 
Z6 = 0 

DO 80 I = 1,11 
Z1 = E(I) 

Z2 = F© 

Z3 = Z5 
Z4 = Z6 

Z5 = Z1*Z3-Z2*Z4 
Z6 = Z1*Z4+Z2*Z3 
80 CONTINUE 

Z1 = Z5 
Z2 = Z6 
T = T4*R9(N) 

Z3 = COS(T) 

Z4 = -SIN(T) 

Z5 = Z1*Z3-Z2*Z4 
Z6 = Z1*Z4+Z2*Z3 
PI = P1+S2*Z5/R9(N) 

P2 = P2+S2*Z6/R9(N) 

40 CONTINUE 

PZ(K) = SQRT(P1*P1+P2*P2) 
WRITE(NOU,3) PZ(K) 

3 FORMAT(3X,’PZ = ’,F9.5,/) 

B=B*T6 

G=Gi>T6 

D=Drr6 

WRrrE(NOO,4) D,PZ(K) 

4 FORMAT(3X,F4.1,3X,F9.5) 

20 CONTINUE 

END 
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